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ABSTRACT 


Experiments were conducted to determine if free-stream 
turbulence scale affects separation of turbulent boundary layers. 
In consideration of possible interrelation between scale and 
intensity of turbulence, the latter characteristic also was 
varied and its role was evaluated. Flow over a 2-dimensional 
airfoil in a subsonic wind tunnel was studied with the aid of 
hot-wire anemometry, liquid-film flow visualization, a Preston 
tube, and static pressure measurements. Profiles of velocity, 
relative turbulence intensity, and integral scale in the boundary 
layer were measured. Detachment boundary was determined for 
various angles of attack and free-stream turbulence. 

The free-stream turbulence intensity and scale were found to 
spread into the entire turbulent boundary layer, but the effect 
decreased as the airfoil surface was approached. When the 
changes in stream turbulence were such that the boundary layer 
velocity profiles were unchanged, detachment location was not 
significantly affected by the variations of intensity and scale. 
Pressure distribution remained the key factor in determining 
detachment location. 
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INTRODUCTION 


This investigation was undertaken to learn if the scale of 
turbulence in the free stream has any influence upon the location 
at which boundary layer separation occurs in an adverse pressure 
gradient. Both free-stream turbulence intensity and scale 
recently have been reported to influence local skin friction 
coefficient (refs. 1-7). Thus, in view of the linkage of that 
factor with separation location provided by Townsend's analysis 
(ref. 8), it is implied that both free-stream turbulence 
intensity and scale may affect separation. If so, there would be 
concern about the varying turbulence scales in wind tunnels, air- 
breathing propulsion units, and numerous other internal flows. 
In addition, the possibility of optimizing flow-manipulating 
devices such as vortex generators would seem feasible. Flow 
separation and related aerodynamic phenomena are of great 
importance in many fields of fluids engineering. 

In view of the prevalence of turbulent boundary layers in 
current practical situations, and also in consideration of the 
expected need to investigate turbulence intensities that would be 
rather high in the context of wind tunnel flows, this 
investigation has been focused largely on turbulent boundary 
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layers. Although laminar separation can occur at low Reynolds 
numbers, in circumstances of practical interest it will 
frequently be followed by transition of the unstable separated 
shear layer and subsequent reattachment as a turbulent boundary 
layer. The latter may separate still further downstream if it is 
exposed to a sufficiently strong positive pressure gradient. 
Thus, it is possible to have both laminar and turbulent 
separations within rather short streamwise distances. 

This investigation was confined to subsonic, essentially 
two-dimensional flows. The evidence for a free-stream turbulence 
scale effect on skin friction seemed somewhat inconsistent (cf. 
refs. 3 and 4) at the time this work was planned, and it was 
thought that the simplest conditions of practical relevance 
should be chosen for study. The approach was, first, to 
determine if free-stream intensity and scale variations affect 
detachment and then to define the parameters that controlled the 
separation process. A si±>sonic wind tunnel with its associated 
pressure transducer and hot-wire anemometry systems and a 2- 
dimensional airfoil comprised the principal experimental tools 
used in this investigation. 

It has been necessary to conduct an extensive experimental 
program to measure free-stream turbulence characteristics prior 
to investigating the influence of those parameters on separation. 
This report includes a summary of those data, which have been 
more fully covered in refs. 9 and 10. Also in connection with 
this research, there was an opportunity to acquire data on 
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laminar separation bubbles, i.e., laminar boundary layer 
separation followed by transition and reattachment. References 11 
and 12 stemmed from that work. 
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TECHNICAL BACKGROUND 


The Separation Zone 

Following Simpson (ref. 13) and others, separation is here 
taken to mean the entire process of breakaway and breakdown of a 
boundary layer from its bounding surface. The mean separated 
flow may be two- or three-dimensional. Although vanishing wall 
shear stress and flow reversal are usually associated with 
separation, these are only typical of two-dimensional separation 
with steady free-stream flow, and it is incorrect to assume that 
they are characteristic of three-dimensional separations or 
unsteady free streams. 

It is essential to note that separation, like boiandary layer 
transition, occupies a spatial zone. This can have important 
implications in connection with the interpretation of patterns in 
liquid surface films or boundary layer profiles measured with 
probes in a separation zone. "Separation" usually is presented 
without qualification as a single-point measurement. The 
following sketch indicates the nature of two-dimensional 
separation with steady free-stream flow, which is the category 
discussed in this report. 
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The Separation Zone 
(Reference 13) 



ID = incipient detachment (1% instantaneous backflow) 

ITD = intermittent transitory detachment (20% instanteineous 
backflow) 

TD = transitory detachment (50% instantaneous backflow) 

D = detachment (mean shear stress = 0) 


One would expect that lengths between the defined points may 
vary with the conditions in a given case. However, there does 
not seem to be much information on this particular aspect of the 
problem. 

The lack of data from probing within separated flow regions 
is repeatedly remarked in survey papers, c.f. ref. 13. High 
turbulence and reversals in mean flow direction, which are 
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characteristic of the separation zone, make the hot-wire 
anemometer unreliable despite its rapid response. Pitot tubes 
also suffer in such environments, and the lag in response of 
their transducers is an added problem. Double-headed pitot 
probes have been used in separation zones to gain directional 
capability, but their other problems, including obstruction of 
the flow, remain. Much can be learned from flow visualization 
techniques, but extremely complex flow fields often are generated 
in separation zones, and the interpretation of surface liquid 
film or "paint" patterns and other visual evidence is not always 
certain . 

Turbulence Intensity 

Free-stream relative turbulence intensity throughout this 
discussion is taken to be 

t 4 = /?/U„ = u 7U<„ (1) 

where 

u^ = mean square of the fluctuating component of the mean 
free-stream velocity, 

A cartesian coordinate system with x in the streamwise direction 
and y normal to both x and the wall is adopted. Orthogonal to 
these is a z-axis considered to lie in the spanwise direction. 
As discussed in ref. 14, a hot wire normal to the x-axis is 
exposed to fluctuating velocity components along all three eixes, 
but under the conditions prevailing in an experiment such as the 
present one, the wire output is dominated by the velocity along 
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the x-axis. 


Turbulence relative intensity within the boundary layer is 
defined as 

T' = u7U (2) 

where 

U = mean streamwise velocity at height y within boundary 
layer 

Turbulence Scales 

To describe turbulent motion quantitatively, in addition to 
relative intensity, the concept of turbulence scales has been 
formulated, (cf. refs. 15-16). Scales in time and in space, as 
well as scales in longitudinal and transverse directions, are 
considered. 

Defining the scales of turbulence usually involves 
correlations between velocity components separated in either time 
or space. In either case, the correlation will be a function of 
time lapse or distance between the points. Distinctions are made 
between small fluid elements (fluid "particles") and larger 
elements (fluid "lumps"). The former are comparable to the micro 
scales, the latter comparable to the integral scales. Eulerian 
and Lagrangian descriptions are useful. In the first, the 
variation of a property with respect to a fixed coordinate system 
is considered; in the latter, the variation of the property 
connected with a given fluid particle or fluid lump is 
determined . 
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Traditionally, the integral scale usually has been used 
in efforts to determine the role of turbulence scale in 
connection with drag, heat transfer, separation, transition, etc. 
However, other scales in the spectrum of turbulence have been 
identified, and several recent reports of studies of the 
influence of free-stream turbulence on boundary layers present 
discussions based exclusively on a dissipation length scale 
parameter (e.g. refs. 5-7). The latter is attractive in that it 
is generally easier to measure. 

With regard to the scale of most influence on boundary layer 
separation, the issue would appear to be whether gradient-type 
transport of turbulence properties, related to the micro- or 
dissipation lengths, or transport by bulk movement related to the 
integral scales is the dominant process. One is attracted to the 
idea that bulk movement, and therefore integral scale, is the 
most likely candidate. The almost universal choice of early 
investigators of scale effect on boundary layers has been the 
integral scales (e.g., ref. 17). Numerous other examples of this 
choice could be cited. However, in 1974, Bradshaw (ref. 18) 
proposed that a dissipation length scale parameter be used for 
interpretation of the effects of turbulence scale upon boundary 
layers, and this is the path followed by several of the recent 
investigators. According to ref. 18, the longitudinal 
dissipation length scale is proportional to and of the order of 
the integral scale. Both scales were measured in the present 
case, and each is defined here. 
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Integral Scale; 


References 15-16, for example, contain detailed discussions 
of this topic, so only the particular approach followed in the 
present work is described. If the fluctuating velocity component 
u is measured at times t and (t + r), a correlation coefficient 
may be defined as 

Ry = u(t) • u(t + r)/[u^(t) • u^(t + r)]°»® (3) 

When is determined for a range of time delays, a curve such as 
shown in Fig. 1 is produced. A quantity defined as 


t = R^ dt (4) 

is multiplied by the mean velocity, U, to give a length 1^ which 
represents the typical x dimension of the energy-containing 
eddies in the flow. This is the integral scale used hereafter in 
this report. Considerations of other coordinate directions and 
velocity components leads to the other integral scales. In 
practice the integration in Eq. (4) is terminated where R^ first 
becomes zero. 


Dissipation Length Scale Parameters ; 

A streamwise dissipation length scale parameter is defined 
following refs. 15 and 18, viz.. 

Lx = - {'^) ^ • V[U<„( d^/dx) ] ( 5 ) 

The longitudinal turbulence energy produced by grids typically 
decays according to a relation such as 

= C(X-Xo)/M (6) 

where u^ = mean u^ value 
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C = a constant 


M = grid mesh size 

X = distance downstream from grid and 

Xq = a constant 

There is frequent need for an effective starting length to be 
used in order to fit Eq. (6) to experimental data, thus X^ is 
present to accommodate that need. When Eq. (6) is differentiated 
and combined with Eq. (5), it is found that 

Lj^/M = nC“®«®/^[(X-XQ)/M]^"°*®/^ (7) 

It is interesting to note that this parameter will increase with 
X/M if n > 0.5 but decrease if n < 0.5. One supposes that the 
value of n is affected if turbulence originates downstream of the 
grid. For example, turbulence may be generated in a diffuser or 
at a wind tunnel fan and propagate upstream in subsonic flow. 
When that happens, the rate of turbulence decay is altered. The 
possibility of turbulence "feed-in" downstream of a grid cannot 
be ignored for many subsonic wind tunnels. 

The boundary layer on an airfoil in a wind tunnel was 
studied in this research. Therefore, free-stream turbulence 
characteristics were measured at a point on the tunnel centerline 
immediately upstream of the station where the airfoil leading 
edge was later located. A catalog of T^^, Ijjoof values was 

compiled for five turbulence-producing grids located at various 
distances upstream of the test section. Examples of these data 
are shown later. 

The previous investigations of the combined effect of both 
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T' and or on local skin friction or Stanton number (cf. 
refs. 1-7) led to the several empirically determined correlation 
parameters listed below. These parameters collapsed or 

AS^/S^q data to single curves with moderate scatter, within the 
scope or range of conditions covered in the referenced studies: 
100 (uVUe)/[(I^/ 6 ) +2] ... ref. 5 

100 ( uVUg)/[ (Ljj/ 6 ) + 2] • [3 exp(-Re 0 / 4 OO) + l] ... ref. 6 
When Castro (ref. 7) extended the AC^/C^q measurements to lower 
Re 0 than previously covered he achieved the correlation of his 
and Hancocks' (ref. 1) data by plotting 
(ACf/Cf^^) { 1 + 10 /[f 2 exp(-Re 0 / 4 OO)] } 
as a function of F, where 

F = 100 (uVUg) / [(0.5Lj^/6) + 2.5] 

In all cases 6 is defined to correspond with U/Ug = 0.995. 

Although the above parameters vary, they all suggest that 
the effects of T' and or are opposite, i.e., increases 
with T' but decreases as or increase. 

A few words are necessary regarding the relationship of 
and L„. In ref. 18 it is remarked that 
Lx “ 4.5 Ijj 

while, in ref. 6 , it is noted that 

Lx = 1*5 Ix- 

In the experiments to be described, neither of the above 
relations holds. In the case of ref. 18, it seems that reference 
is made to grids installed upstream of the working section in a 
duct of constant area. In ref. 6 , it is clear that the grids 
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were mounted upstream of a contraction in area preceding the test 
section. Both configurations were used in the present work. 

Interaction of Free-Stream and Boundary Layer Turbulence 

It is perhaps premature to present this topic, because it 
will have to be taken up again when results of this investigation 
are discussed. However, for background to the experimental data 
as they are presented, some of the physical concepts applicable 
to the present study should be mentioned. These are drawn 
principally from refs. 4, 7, 13, and 18. 

Meier (ref. 4) comments that "Maximum values of the skin 
friction coefficient at T« = const, were obtained at length 
scales in the order of the boundary layer thickness." Simpson 
(ref. 13) observes that "the backflow region is strongly 
dominated by turbulent fluctuations that are greater than, or at 
least comparable to, the mean velocities . . . ." Further, he 
states that there is no location with backflow all of the time. 
This is significant in connection with efforts to measure 
velocity profiles in separation zones. 

In the introductory part of ref. 18, Bradshaw suggests 
tentatively, based on various evidence, that "free-stream 
turbulence reduces the ability of the existing shear layer 
turbulence to transfer momentum or scalars, and/or that the 
shear-layer turbulence reduces the mixing ability of the free- 
stream turbulence that adjoins it." At another point, he asserts 
that "it appears that large-scale free-stream turbulence has no 
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direct effect on turbulence structure. . . Also, "If the 
intensity and length scale of the free-stream turbulence are of 
the same order as those of the shear layer - as they must be to 
produce significant effects. . . Lastly, "It is demanded by 
logic, and confirmed by experiment, that the effect of free- 
stream turbulence is first felt by the outer layer." 

Castro (ref. 7) writes that, "It appears that the addition 
of free-stream turbulence significantly reduces the Reynolds 
number effects on u^ and there is even a suggestion that at the 
highest value of FSTP (free-stream turbulence parameter) . . . 
the effect is reversed." The FSTP referred to is 

(u7Ue)/[(Lj^/6) + 2] 

This section has been included to provide background to the 
remainder of the report. The linkage between the earlier work, 
most of which reported the effect of free-stream turbulence 
on Cf or S^, and separation is offered by Townsend's analysis 
(ref. 15). It is shown there that is a factor affecting 
separation location. That relationship is made use of later when 
the present data are discussed. 
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EQUIPMENT and TECHNIQUES 


Wind Tunnel 

The Vanderbilt University Engineering School Wind Tunnel is 
shown in Figures 2 and 3. A top view of the tunnel is not shown 
because the structure is of square cross section upstream of the 
diffuser, and the latter component varies from square to circular 
as it extends downstream. Except for the fan section, the tunnel 
shell is made of fiberglass, giving a very smooth interior sur- 
face. A summary of flow conditions is presented in Table 1. 

The smaller, higher speed, 0.4-m section is the normal test 
section, and all subsequent comments will pertain to that area 
unless specific reference is made to the 1-m section. 

A shroud, honeycomb, and two screens are installed at the 
wind tunnel entrance. The honeycomb is fabricated of plastic 
impregnated paper having a thickness of 0.2 mm. There are two 
layers making up the honeycomb. In each layer the cells, or 
openings, are approximately triangular in cross section, with a 
base of 9.5 mm and height of 8 mm. Streamwise dimension of the 
assembly is 64 mm, giving a length-to-hydraulic diameter ratio of 
approximately 7. 

Two fiberglass screens with 6.3 mesh/cm (16/in.) are 
installed to reduce turbulence. Additional screens may be 
inserted if it is desired to study phenomena requiring very low 
free-stream turbulence. Relative turbulence intensity in the 
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test section is approximately 0.16 percent with the present set 
of screens. Figure 4 is a plot of T« as a function of free- 
stream velocity. 

Velocity distribution laterally across the 0.4-m test 
section, measured by hot-wire and pitot-static probes, is 
constant to within ± 1 percent in the central 0.3 m (12 in.) at a 
mean velocity of 60 m/s (200 fps). Longitudinal velocity 
gradient in the test section is virtually zero at all of the 
usual test conditions. 

A flexible coupling at the downstream end of the diffuser 
reduces vibration that would be caused by the centrifugal fan and 
motor. The laboratory floor is of thick concrete and quite 
rigid. However, at higher speeds there is a small amount of 
vibration in the tunnel shell, so sensitive probes are mounted 
onto the laboratory floor or to heavy structural steel beams 
above the wind tunnel. Most of the flow surveys discussed later 
were made with hot-wire or pressure probes mounted to the 
overhead structure as shown in Fig. 5 

Figure 6 shows the airfoil with the spoiler used in this 
investigation, and Fig. 7 shows the fiber-optic, surface-sensing 
device that enables the delicate hot-wire sensor to be positioned 
within a quarter of a millimeter or less distance from a solid 
surface without being broken by contact. The latter is a fiber- 
optic, light-emitting auxiliary probe (A) which is mounted on the 
same micrometer system that positions the flow probe (B). Light 
from a photo-flood light outside of the tunnel passes through a 
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fiber-optic filament and is emitted from the end of the tube (A), 
which is adjusted so that it is perpendicular to the model 
surface (C). When the end of this tubular probe contacts the 
model surface, the circular spot of light beamed from its open 
end (D) is extinguished. This provides a sufficiently precise 
and repeatable indication of hot-wire distance from the model 
surface when the end of tube (A) touches the surface and the 
displacement between hot-wire flow probe and light-emitting probe 
tip has been determined. The latter dimension is measured in the 
present case by a telescope equipped with a micrometer and 
movable cross hairs. In practice, the two coupled probes are 
displaced so that the very fragile hot wire sensor does not 
contact the model surface when the light probe touches that 
surface. The principal advantage of this technique for indexing 
probe position is that it is usable when the model surface is not 
electrically conducting and an electrical signal of surface-to- 
probe contact cannot be obtained. 

Pressure and Hot-Wir e/Hot-Fi Im Systems 

The Vanderbilt Wind Tunnel is equipped with a Scanivalve 
Corporation DSS 24C MK3 electronic switching-valve pressure 
transducer system. When this is used in conjunction with a 
manually-operated Norgren Fluidics 10-port switch, up to 34 
pressure measurement channels are available. Because the 
Scanivalve pressure transducer range is too high for accurate 
results in the present program, a MKS Baratron Model 170M-6C 
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precision pressure transducer with a full-scale range of 100 mm 
Hg was used. 

Three channels of hot-wire instrumentation are available. 
These include one Thermo Systems, Inc. (TSI) Model 1050 and two 
TSI Model 1010 anemometer systems. Each anemometer channel has a 
linearizer, a TSI Model 1052A being used with the 1050 anemometer 
and TSI Model 1072 linearizers with the 1010 anemometers. A TSI 
1015C Correlator is available for correlation of signals from the 
anemometers. However, it was more convenient and accurate to 
utilize digital signal processing techniques made possible by 
availability of an analog-to-digital convertor and an Apple lie 
computer connected to VAX and DEC mainframes. 

The hot-wire voltage measurements, D.C., rms, and mean 
square, were obtained with a TSI Model 1076 digital voltmeter. 
This instrument may be operated with time constants of 0.1, 1.0, 
10 and 100 seconds. 

Data acquisition was accomplished with an Interactive 
Structures, Inc., Model AI13 analog-to-digital converter coupled 
to an Apple lie computer. Experiments were conducted to 
determine the effects of sampling time and frequency, and the 
optimum schedule was selected. Computer programs written to 
process the hot-wire data and compute turbulence properties are 
in Appendix I of this report. 

The hot-wire/hot-film probes were of convention design. The 
hot-wire sensors were platinum coated tungsten on gold plated 
needle-type supports. Diameter of sensing element is 3.8 y 
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(0.00015 in.) and the length is 1.27 mm (0.050 in.). 

Corresponding dimensions of the hot-film probe sensors were 51 y 
(0.002 in.) and 1.0 mm (0.040 in.). The sensor films were 
platinum with an alumina coating. Fused quartz substrates and 
gold plated needle-type supports were used. A photograph of a 
hot-wire probe may be seen in Fig. 7. 

Although some early measurements were made with hot-film 
probes, it was discovered that those gave erroneous (high) u' at 
free-stream speeds above roughly 40-45 m/s (140 fps). This 
problem was not studied in detail, once the trouble was isolated 
to the probe type. However, it is apparent that the much greater 
diameter of the hot-film sensor led to much lower frequencies of 
vortex shedding (cf. ref. 19), and it seems likely that this may 
have affected the readings. Probe vibration was also suspected, 
but the addition of stiffening fillets to the sensor-supporting 
needles did not eliminate the trouble. It is noteworthy that a 
similar failing was found later when using a "boundary-layer" 
type of probe with long, curved, needle— like supports for the 
sensor, even though these had a hot-wire sensor. In that case, 
it seems most probable that the long support needles must have 
vibrated. Changing to stiffer type of probe seen in Fig. 7 
apparently avoided the problem. All hot-wire data presented in 
this report were obtained with probe configurations that 
performed best and seemed free of spurious influences after a 
period of trial testing. Experiments were performed to confirm 
that the support strut of the probe shown in Fig. 7 did not 
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affect the sensor signal. 

It was necessary to make a small correction to the hot-wire 
voltages when the wire was very close (< 1 mm) to the airfoil 
surface. Heat transfer to the cooler airfoil caused hot-wire 
voltages to falsely indicate higher velocities. An example of 
this effect is shown in Fig. 8. The correction procedure 
followed ref. 20. 

Another type of error becomes important when turbulence 
intensity levels are above the order of 10%. Near the airfoil 
surface in the separation zone, turbulence intensities of 10 to 
50% were recorded. Following Bradshaw (ref. 14) a maximum error 
in T' of approximately +11% is estimated for the most adverse 
conditions for which data are given in this report. Inasmuch as 
that one case is only used as a qualitative example of high T' 
encountered in the separation zone, and all other corrections 
would be within the uncertainty band of the data, ' no corrections 
for high -turbulence error were made. 

Skin friction measurements were made with the Preston tube 
shown in Fig. 9. Outside diameter at the entrance was 7.1 mm 
(0.028 in.). The pressure was read with the MKS Baratron system. 
The procedure followed ref. 21, and it was found that the 
parameters in the present case fell within the range covered in 
that reference. Effects of pressure gradient on Preston tube 
data were evaluated following ref. 22 and foimd to be negligible 
for the measurements reported here. 
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Detection of Detachment By Liquid Films 

Liquid films of various solutions were painted onto an 
airfoil and tested in the tunnel. It was found that a mixture of 
the following approximate proportions (by volume) served our 
purpose ; 

50 percent water [Optimum proportions 

35 percent dry powdered poster paint varied with tunnel 
10 percent glycerine speed. This was 

5 percent Kodak Photo Flo. typical.] 

This solution was brushed onto an airfoil surface that had been 
wetted with a film of soapy water. Air flow was immediately 
started, and the drying liquid film was observed until an 
unchanging indication of "detachment" was seen. 

Photography of the dried paint had to be done with the aid 
of a pair of plane mirrors because of the confined area and 
difficulty of access created by other apparatus. However, this 
technique was quite successful and only involved the extra work 
of installing and removing the mirrors between wind tunnel runs. 

Photographs and tabulated data are presented later where the 
interpretation of the paint results is addressed. In general, 
very satisfactory information was obtained, both regarding 
separation and the quality of two-dimensionality of the flow over 
the airfoil. The streamwise streak lines visible in the dried 
paint are seen to be parallel to the free stream and normal to 
the spoiler on the airfoil. 


20 



The Airfoil 


An airfoil was chosen to be the body on which to study 
separation because it is obvious that one of the principal areas 
of application of this research is that of airfoil and wing 
design. Figure 10 is a sketch of the airfoil. When installed, 
it nearly spanned the test section of the wind tunnel. End 
plates of relatively large size were fitted initially, but the 
end-plate area above the airfoil surface was soon cut away to 
avoid boundary layer development in the corners between the 
airfoil upper or working surface and the plate. Leaving end 
plates extending well beyond the lower surface seemed sufficient 
to prevent excessive air spillage into the region of the boundary 
layer surveys on the upper surface, and flow-visualization by the 
liquid-film or paint technique revealed a larger span of two- 
dimensional flow without the interference of the boundary layer 
on the upper area of the end plate. 

Orifices for pressure measurement were located at the 
nominal stagnation point and every 1.27 cm (0.50 in.) downstream 
on the upper surface. Angle of attack could be varied, but most 
experiments were conducted with an angle of 1 deg. It was found 
that better, i.e., two-dimensional upper surface flow existed if 
a relatively low angle of attack was maintained and separation 
was caused by a spoiler, as shown in Fig. 10. This also imposed 
less aerodynamic load on the airfoil and mountings, so that there 
was less movement under load and more precise probe positioning 
was possible. 
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Grids and Effects of Wind Tunnel Contraction Section 


Five different grids were used in the course of this 
investigation. Figure 11 gives the dimensions. It will be noted 
that Grids 1-3 were mounted in the tunnel upstream of the 
contraction that separates the 1-m and the 0.4-m test sections. 
Grids 4 and 5 were mounted at the downstream end of that 
contraction, i.e., at the upstream end of the 0.4-m test section. 

Streamwise positions of Grids 1-3 were variable. Continuous 
tracks and movable mounting clips made that possible. However, 
Grids 4 and 5 were located at only one streamwise station as 
already specified. Figure 12 shows some of the grids. 

It is shown in refs. 9, 10, and 23 that the effect of a 
contraction in cross section area upon stream turbulence is 
dependent on the character of the turbulence entering the 
contraction. Passage through a contraction between grid and hot- 
wire probe stations also alters the relationship between 1^^ and 

In the present case, the wind tunnel contraction area ratio 
is 6.2 (geometric) or 6.67 (based on measured center line 
velocities). Integral scales were measured at the entrance and 
at the exit of the contraction section, in a brief investigation 
of this phenomenon, and it was found that a magnification of 
integral scale occurred at all velocities and turbulence 
intensities. Figure 13 presents the measurements of Ix2/^xl ^ 
function of T^. I^^^ does not appear to be a factor within the 
range of these data. Two empirical relationships roughly fitting 
the data for 1.5 < Tj'% < 6 are 
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or 


( 8 ) 

(9) 


Ix2/Ixl ^ 2.67/T{% 

W^xl "" 3.71/(T{%)“*" 

Figure 13 also shows these equations but it has to be emphasized 
that they are only based on these limited observations. 

It may be significant that the contraction section where 
these data were obtained has a streamwise length that is 0.775 
times the lateral dimension of its entrance and is of square 
cross section. The change in turbulence intensity as flow passes 
through an area contraction has been found to be largely 
dependent on the area or velocity ratio alone (cf. ref. 24 ), 
although more complex relationships have been published (ref. 
25). However, the possibility of the rate of change of velocity 
being a factor in determining the ratio of Ix2/^xl should not be 
overlooked. 

Turbulence intensity ratios (T^/T^) corresponding to the 
data in Fig. 13 average 1/6, which is very nearly identical to 
the wind tunnel contraction ratio, A 2 /Aj^. Figure 14 demonstrates 
that this result is consistent with data for some NASA tunnels 
given in ref. 24. On this subject, it should be noted that some 
decline in T' would occur over a length of constant area duct, so 
that all of the decrease in T' shown in Fig. 14 cannot be 
credited to the contraction in area alone. 

Table 2 lists the T<^ and 1^00 values that made up the menu of 
free-stream turbulence conditions available for selection in the 
investigation of separation. Having a range of T^^ with 
approximately equal to turbulent boundary layer thickness on the 
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downstream one-half of the airfoil chord was a goal. The damping 
effect on exerted by the contraction section made it difficult 
to obtain higher T(j^ using grids 1-3. Therefore, grids 4 and 5 
were fabricated and installed downstream of the contraction. 
This produced higher levels of T« but the position 30 grid mesh 
lengths upstream of the airfoil leading edge for grid 4 and 15 
mesh lengths for grid 5 was viewed with some apprehension because 
of possible flow nonuniformity. This concern was abated when a 
survey across the stream in the test section showed very 
acceptable uniformity of turbulence intensity. These data are 
presented in Fig. 15. The separation data obtained when grids 4 
and 5 were used appear to be consistent with the results obtained 
when grids 1,2, and 3 were used. 
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EXPERIMENTAL RESULTS AND ANALYSIS 


Pressure Distribution and Transition 

Figure 10 shows the airfoil and spoiler arrangement. The 
pressure distribution on the airfoil for angles of attack of 1, 
3, 5, and 9 deg and free-stream speed of 5 0.5 m/s are presented 
in Figs. 16 a-d. 

The location of boundary layer transition with no trips was 
determined by moving a hot-wire probe longitudinally along the 
surface at y = 0.15 mm and recording T". The results are given 
in Fig. 17 for two unit Reynolds numbers. With Re„/m = 3.15 x 
10® the values of ROoot 6.07 x 10® for the beginning of 
transition and roughly 9.6 x 10® at the end of transition. When 
Rego/m = 4.00 x 10®, transition began at = 7.12 x 10® and was 
completed at Re^o = 8.74 x 10®. Because it was preferable to 
conduct most of the experiments with U<,o = 50.5 m/s, which 
corresponds to Re„/m = 3.15 x 10®, trips were used to assure a 
fully turbulent bovindary layer at the beginning of the adverse or 
positive pressure gradient on the airfoil. Again, a hot-wire 
probe was used to verify the desired fully turbulent boundary 
layer while using a trip of minimum height. 

Detachment Locations 

Figures 18 and 19 are examples of the data on separation as 
determined by the liquid-film technique. The small scale shown 
on the right-hand side in all overhead photos is pushed against 
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the upstream face of the spoiler. Numerous repetitions of the 
liquid-film measurements demonstrated excellent consistency of 
results. Before proceeding to discuss these data it is necessary 
to more precisely relate these observations to the terms in the 
sketch of the separation zone given in the Introduction. These 
photographs and numerous other examples of the liquid-film all 
display the major features sketched below. 



(d)- 


spoiler 

(a) = region approximately 1 cm in streamwise extent 

believed to be an area of intermittent backflow 
containing the ID and ITD boundaries defined earlier, 
following ref. 13. Paint in this area apparently is 
swept downstream where it contributes to formation of 
(b). 

(b) = highly repeatable, narrow, heavy band of paint 


26 


believed to indicate near 50% instantaneous backflow, 
i.e., the TD station of ref. 13. Its prominence in the 
photos is partly due to its light-reflecting property. 
Detachment appears to occur essentially at the 
downstream edge of this band, making the TD and D 
stations very close together in this case. 

(c) = an area several cm in streamwise extent where 

relatively low and random surface shear is evident. 
When the paint in this region was not overflowed by 
paint from other regions, it remained smooth and uni- 
form while it dried. 

(d) = highly repeatable, roughly 1-cm region at the foot of 

the spoiler where the surface is cleansed of paint by 
energetic vortical flow. The paint is swept upstream 
to the rear boundary of region ( c ) . 

Kot-wire and small impact-pressure probes were used to 
survey the regions (a-d), but the large magnitudes of fluctuating 
velocities and intermittency of velocity directions make probes 
of these types ineffective for determining details of flow near 
the surface in separation regions. (Use of a special laser- 
Doppler anemometer for this type of measurement is reported in 
ref. 26.) The typical, orderly velocity profiles with reversed 
flow near the solid surface that are often sketched when 
separation is discussed qualitatively do not adequately portray 
the intermittency and transitory patterns of real separation. 

Further in regard to interpretation of the mean detachment 
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station, typical surface static pressure distributions across the 
separation zone are shown in Figs. 20 a-b. The "detachment" 
label in these figures is placed where the paint band (b) was 
located. It will be noted that the pressure coefficient did not 
reach a constant value at the designated detachment, but it did 
become constant a short distance downstream. This feature is 
common to many similar examples of pressvire distributions with 
separation present. Thus, it seems justified to conclude that 
the dark paint band (b) in the photographs marks the upstream 
boundary of detached flow. 

The Influence of Ti and 

Table 3 and Figs. 21-22 present detachment station, at 
corresponding T® and l^oor ss determined in this investigation of 
subsonic, turbulent boundary layers in typical adverse pressure 
gradients on the suction side of a two-dimensional airfoil. 
Angles of attack and free-stream velocities are indicated with 
the data in Table 3. Figures 21 and 22 correspond to a = 5 deg 
and 50.5 m/s. Tables 4 and 5 present data selected so that 
either T<^o or nearly constant while the other parameter 
varied. This format makes it easier to see how or if each 
turbulence parameter independently influenced detachment point. 
The airfoil was at a = 5 deg. and 50.5 m/s conditions when the 
data of Tables 4 and 5 were recorded. 

The range of T^ adequately represents wind tunnel flows and 
the higher turbulence of some other types of internal flows. As 
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planned, the free-stream integral scales are approximately equal 
to the boundary layer thicknesses upstream of the separation 
zone. It will be recalled that several previous investigators 
have concluded that this should lead to the maximum influence on 
surface shear stress. As Table 3 shows, the boundary layer 
thickness grew rather rapidly between the mid-chord and 
detachment stations. Thus, generally fell between the values 
of 6 at X = 22.1 cm and x = 29.7 cm in these experiments. 
The ratio given in the tables is based upon §^995 at x = 
29.7 cm, which is approximately the detachment station in all 
cases. If 6 at x = 22.1 cm were used, all ratios would be 
approximately doubled because the boundary layer thickness grew 
rapidly toward detachment. 

The evidence in the tables and figures seems decisive; in 
the range of T« and 1 ^^^ investigated, neither parameter had a 
significant effect on x^j. All measurements were repeatable to a 
very satisfactory degree, and no anomalies were observed. When 
confronted with this result, the decision was made to "back- 
track" and briefly check the influence of free-stream turbulence 
on skin friction. 

Measurements of 

Earlier investigators (e.g., refs. 2 and 18) have reported 
on this, and it was seen as a way to determine if the present 
experiments were consistent with the related previous 
investigations. Failure in this test would point to some defect 
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in the present experiments. Therefore, a Preston tube was used 
to determine Cj at mid chord for several values of . Results 
are presented in Fig. 23. 

Inspection of figures presented in refs. 6 and 7 reveals 
that the data on = f(T^ ) typically scatter across the 

band drawn in Fig. 23. It has already been remarked in the 
previous section that there is agreement among refs. 3-7 that 
either integral length scale or dissipation length parameter also 
is a factor that should be incorporated into the correlation of 
Cj/Cjq with free-stream turbulence. However, the best form for 
this correlation parameter seems to be uncertain, and in any 
event, it apparently will be a weaker factor than . 

Therefore, Fig. 23 implies that the present experiments are 
consistent with earlier results, and no suspicious discrepancy is 
evident . 

Comparisons With Prediction Methods 

The final step in confirming the detachment results was the 
comparison of the measured with predictions by Stratford's 
(ref. 27) method. That result is seen in Table 6. When making 
the calculations of x^^, the experimental Cp distributions were 
used. 

The designations Strat. 1 and Strat. 2 refer, respectively, 
to the version of Stratford's method originally presented in ref. 
27 and the modified method proposed by Cebeci, et al. in ref. 28. 
The modified version incorporates a small change in the critical 


30 



value of Stratford's separation criterion. In both cases, the 
initial length of laminar flow and favorable pressure gradient 
that existed on the airfoil in these experiments is taken into 
account. When the modification recommended by Cebeci, et al. is 
made, Table 6 shows excellent agreement between Strat. 2 and 
experiment. 

Inasmuch as the surface shear stress or Cf appears directly 
in Townsend's method for predicting (ref. 8), a brief 
examination of the effect of increased Cf caused by elevated 
free-stream turbulence was carried out. This was based on the 
faired curve for the "present data" in Fig. 23. That curve lies 
close to the Simonich and Bradshaw curve within the range of T^ 
in Fig. 23. 

The calculations by Townsend's method generally agreed 
closely with the results obtained by Stratford's method, as 
demonstrated in Table 7. However, in a number of cases 
Townsend's procedure failed to predict detachment, and 
Stratford's method gave better results. Use of the modified 
Townsend method proposed by Hahn, et al. (ref. 29) was even less 
successful in those cases. When either or corresponding 
to elevated turbulence (from Fig. 23) were used in the Townsend 
procedure for predicting x^j, the difference was insignificant. 

Once again, it should be remarked that the experimental 
pressure distributions were used for all of the x^ calculations. 
The Townsend method actually forecasts the Cp at detachment. 
When detachment occurs at a Cp lower than predicted, the 
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predicted value is not attained and there is simply no solution. 


Interaction of Free-Stream and Boundary Layer Turbulence 

In consideration of the basic data, the overall consistency 
of all results, and the generally very good agreement found where 
there were check points provided by related results from other 
investigators, the earlier statement that there is minimal 
influence of moderate levels of free-stream turbulence on is 
now more strongly advanced. Of the pair of character istics,T.' 
and Ij^cDf neither the usually stronger parameter nor the 

weaker 1^^^ were foxand to have a measurable linkage with at the 
levels investigated. The fundamental factor would seem to be the 
weak penetration of free-stream turbulence into the inner part of 
the turbulent boundary layer. In pursuit of this subject, 
profiles of U/U^, T' , and 1^^ were measured at the stations 
X = 22.1 cm and x = 29.7 cm for different levels of T^ . The 
former station is well upstream of the separation zone, while the 
latter station is in the upstream or early part of that zone. 
The results for x = 22.1 cm are shown in Figs. 24 - 26. 

It will first be noted that change of T« from 6.16% to 
3.16% had no effect on thickness of the boundary layer. The 
effect of that increase in T^ is apparent in Fig. 25, and the 
increase of Ij^^^ also affects the profiles in Fig. 26. Next, the 
profiles of the same quantities at x = 29.7 cm are examined to 
learn how these turbulence properties change as the separation 
zone is entered. 
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Figures 27 - 29 display the profiles of U/U^, T' , and I 
found at X = 29.7 cm. At this station (approximately at 
detachment), the situation is different in some respects from the 
X = 22.1 cm case just discussed. First one notes that changing 
again has an insignificant effect on the velocity profile. 
However, boundary layer thickness has increased markedly in 
comparison to the data for x = 22.1 cm. 

The response of the boundary layer profiles to free-stream 
turbulence intensity and scale shown in Figs. 28 and 29 is 
similar to that found for those two parameters at x = 22.1 cm. 
These figures again show that the feed-in of free-stream 
turbulence characteristics weakens as y decreases and has reduced 
influence near the solid boundary. 

A significant point well illustrated in Figures 25, 28 and 
30 is that T' near the solid surface in the separation zone is 
much greater than the free-stream turbulence intensity. That 
strongly implies that only very high levels of free-stream 
turbulence may have any possibility of affecting separation. Of 
course, very great turbulence will also affect boundary layer 
growth and thereby affect the pressure distribution. The latter 
factor is dominant in regard to detachment and could easily 
overshadow any direct effect of changes in turbulence intensity 
within the boundary layer. 

Figures 26 and 29 show that integral scales within the 
boundary layer had a characteristic profile that contained scales 
lower than free-stream near the wall and greater than free-stream 
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in tile outer part of the boundary layer. There was a consistent 
increase of as increased, but the effect was small at the 
y - values deep within the boundary layer. 

When examining the profiles of some consideration should 
be given to the nature of the autocorrelation procedure and the 
degree of uncertainty. A repeatability of approximately ± 10% is 
indicated by measured on different days under the less 
favorable signal - to - noise conditions corresponding to low 
free-stream turbulence. For a given run or profile in Figs. 26 
and 29, the uncertainty must be less because T' is much greater, 
making signal - to - noise ratio higher, measurements were taken 
within a brief time interval, and only one wind tunnel run was 
involved. However, conditions near the outer edge of turbulent 
boundary layers are such as to impose greater scatter of 1^^ 
values. 

To satisfy the investigators' curiosity, a last effort was 
made to achieve some effect on detachment through free-stream 
turbulence. A horizontal bar of 0.635 cm diam was located 9.5 cm 
upstream of and parallel to the leading edge at a height such 
that the wake of the rod impinged on the airfoil at the highest 
point or "crest" of the airfoil. This did finally cause 
detachment to move appreciably downstream. However, prior to 
accepting this as a direct result of free-stream turbulence, per 
se, the pressure distribution on the airfoil was measured. When 
the pressure distribution was used with Stratford's method, it 
was found that the entire change in x^ was predicted on the basis 
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of the altered pressure distribution alone. Therefore, it can 
not be concluded that the gross increase in turbulence, in 
itself, had any direct effect on x^. Apparently the beneficial 
effects of vortex generators are attributable to their being more 
efficient in redirection of higher-speed air into the inner 
boundary layer than the less-ordered turbulence of the transverse 
rod. In either case, the indirect influence on pressure 
distribution arising from altered boundary layer profiles should 
not be overlooked. 

Prandtl s mixing length concept and the relationship of 
mixing length to viscous shear stress are well known. (See, 
e.g., ref. 30). For 2 -dimensional, parallel flow, this relation 
may be expressed as 

T = |dU/dy I (dU/dy) (10) 

where t = shear stress 
and a = mixing length 

von Kantian 's equation for mixing length is (cf. ref. 30) 

5. = K I (dU/dy)/d^U/dy2| (11) 

It has been found experimentally that H » 0.4. Hence, for typical 
turbulent boundary layers, mixing length is essentially 
proportional to y within the inner part of the bottndary layer. 

Prandtl likens I to the mean free path of kinetic theory of 
gases. In the development of an equation for the mean free path 
for perfectly elastic, spherical molecules of rarefied gases, it 
is found that (ref. 31) 
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X « l/(Na) 

where X = mean free path 


( 12 ) 


N = number of molecules per unit volume 
a = effective cross sectional area of molecules 
Now stretching this obviously simplified model to 
accommodate mixing length and mean free path concepts, one may 
visualize "turbulent lumps" of a diameter proportional to 
integral scale interacting according to the mixing length and 
mean free path concepts. Then, N could be regarded as the 
counterpart of l/I^^ ^ the counterpart of If mixing 

length is analogous to mean free path, it follows from Eqs. (10- 
12) that a and should have similar distributions with y in the 
inner part of the boundary layer. The distributions of Ix^Y^ 
Figs. 26 and 29 seem to support this argument. It has to be 
remembered that local mean velocity U(y) is used with the 
integrated correlation coefficient (Eq. 4) to obtain Ij^, so it is 
not surprising that Ij^ tends to zero as the solid surface is 
approached. Obtaining I^^ in this way probably leads to some 
inaccuracy very near to the wall where T'(y) is relatively large. 

In this connection, values of Ij^ at or just beyond the outer 
edge of the boundary layer do not necessarily equal the free- 
stream values because the local mean velocity may be greater or 
less than the free-stream value. In Fig. 26, = 1.13 U^or and 

in Fig. 29, Ug = 0.96 U„. Thus, the integral scales are seen to 
be approximately = ^x<» (Ug/Uo,) when y ^ 6. 
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CONCLUSIONS 


It has been learned that neither a 20-fold increase in free- 
stream relative turbulence intensity nor a variation of free- 
stream integral scales produced a significant change in boundary 
layer detachment location on the subsonic airfoil studied. The 
integral scales were on the order of boundary layer total 
thickness, and data were obtained in which turbulence intensity 
and integral scale were varied independently. A turbulent 
boundary layer was maintained upstream of the separation zone. 
All intermediate experimental data, such as the influence of 
turbulence intensity upon skin friction coefficient, agree in all 
essentials with results published by others. 

The penetration of free-stream turbulence into the turbulent 
boundary layer was also investigated in an effort to shed light 
on the lack of influence of free-stream turbulence properties on 
boundary layer detachment. Results presented herein clearly show 
how velocity, turbulence intensity, and integral scale profiles 
were affected when the latter two parameters were varied in the 
free stream. In this experiment, the velocity profiles are not 
affected, but the intensity and scale profiles clearly show 
interaction with the free stream. Despite the consistent and 
plausible evidence of the expected penetration of turbulence, 
under the conditions of this research, this interaction did not 
change x^ to any major degree. Such a change in x^ only occurred 
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when a gross increase in free-stream turbulence was brought about 
by placing a bluff obstacle immediately upstream of the airfoil, 
and it appears that the effect on actually should be credited 
to the altered pressure distribution that accompanied the much 
thickened boundary layer. 

It is shown that the influence of free-stream turbulence 
diminished as the airfoil surface was approached. This lessened 
effect near the surface, plus the weak influence of Cj on and 
the dominance of pressure gradient seems to account for the 
insensitivity of to the free-stream turbulence characteristics 
as long as the boundary layer character and thickness 
distribution is unchanged. 

When Stratford's procedure for predicting was used in 
conjunction with the measured pressure distributions and the 
computational modification presented by Cebeci, et al., very good 
agreement with measured x^j resulted. 
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APPENDIX I 


DATA COLLECTION AND REDUCTION PROGRAMS 

The following progreuns were written and used to col- 
lect turbulent velocity data and to calculate the turbulence 
intensity, integral scale, and dissipation length parameter 
from the raw data. The purpose of each program and its 
listing is included. The comments in each program describe 
its operation. 


A. GETSLOWDATA 


This BASIC program for the APPLE lie computer is used 
as indicated in Appendix II to collect the basic velocity 
data needed to calculate the turbulence intensity and in- 
tegral scale. The assembly language subroutine FASTER. OBJ 
is the actual data collection routine. 



QRlGIf^t If 

OF POOR QUAUTt 


10 HOME : CLEAR ! HIMEM: 38144 
20 DIM AX< 15200) .Br.OOO) 

30 Ry. - 15000:MX = 384:MTIME = 54:CT1ME = 10:M1CR0 = 1E04 

40 DM = 6000:BFS = 5000:8107. = 4096:01 = 0 

50 LTAB 5: HTAB 5: PRINT ‘LOADING MACHINE LANGLIA3E SUBROUTINE " 

60 D» « CHR* <4>: PRINT D»i‘BL0AD MONDELAY . OBJ ,A»9500 ‘ 

70 HOME : »JTAB 5: HTAB 10 s PRINT ‘Cl^WIEL flUMBER* ; : HTAB 30: PRINT ‘INPUT 

80 HTAB 16: PRINT *0";: HTAB 31: PRINT •RMS*: HTAB 15; PRINT M5-j; HTAB 

31: PRINT -riKS": V/TAB 11: HTAB 20: PRINT "RANGE*;: HTAB 30: PRINT *C0 
DE" 

90 HTAB 18: PRINT *0 - 5 V ; : HTAB 32: PRINT -0": HTAB 18: PRINT *0 “ I ^ 

HTAB 32: PRINT “1*: HTAB 17; PRINT "0 - .5 V*;: HTAB 32: PRINT *2 
•: HTAB 17: PRINT "0 - .1 V;; HTAB 32: PRINT "3*: HTAB 5: PRINT “ADD 
4 TO GAIN CODE FOR BIPOLAR RANGE- 

100 PRINT : PRINT ; INPUT -INPUT GAIN FOR CH^^EL HO -;G0;C0GX = 0 ♦ 16 * 
G0:C1G7. » 15 + 16 • 61 
no POKE 6,C06X: POKE 7,C167: 

120 HOME : PRINT CHR^ (7): UTAB 15 
130 INPUT -S^PLE TIME? <SEC) -;TIME 

140 INPUT -NUMBER OF SAMPLES? -;NAX: IF NAX < = RX GOTO 160 

150 PRINT -r^IMUM SAMPLES EXCEEDED": PRINT RX;“ SAMPLES TAKEN * :NAX = RX 

160 DR = TIME / NAX * MICRO: IF DR < MTIME GOTO 180 

170 DA = (DR - MTIME) / CT1ME:AH>: = INT (DA / 256) ;ALX = DA - AHX • 256: GOTO 

190 

180 AHX * 0:ALX = 0 

190 SPACE = ((AHX • 256 ♦ ALX) • CTIME + MTIME) / 1000: DR = SPACE • 1000 
195 PRINT : PRINT : PRINT -TIME BETWEEN SAMPLES = *; SPACE ; * msec* 

200 DB * (DM - MTIME) / CTIME;BHX » INT ( DB / 256>:BLX * DB - BKX • 256 

210 POKE 64, AHX: POKE 65, ALX: POKE 66,BHX: POKE 67 , BLX 
220 AX(0) « NAX / 1 28 ♦ 1;BX(0) - 3 

230 PRINT : PRINT : INPUT -READY TO SAMPLE ? -;60^ 

240 POKE 8,1 
250 CALL 38144 
260 BSUM = 0 

270 FOR I = 1 TO MX: BSUM = BSUM + BX(I): NEXT 
280 BAVG = BSUM / M!>; 

290 BVOLT * BAVG « BPS / BI6X 

300 PRINT : PRINT -MKS UOLTAGE == -;BV0LT;- mU":B^)0LT « BVOLT / 1000 

310 PRINT : PRINT : PRINT : INPUT 'SAVE DATA? (ENTER 1 FOR YES, 0 FOR NO) 

- ;60X 

320 IF GOX < 1 GOTO 480 

330 HOME : LnTAB 22: INPUT "ENTER DATA FILE NAME * ; FS 
340 INPUT -DISK DRIVE NUMBER ? " ;AS 

350 PRINT : PRINT : INPUT "DC VOLTAGE = ? (mV) * ;VDC 
355 CALL 38288 

360 CDS « -,D-:NH ^ Ft * CDS ♦ A* 

370 PRINT : HTAB 10: PRINT "SAVING * ;F^ 

380 PRINT DS;-OPEN • ;NS 

390 PRINT DS; -WRITE * ;FS 

400 PRINT 

410 PRINT GO 

420 PRINT DR 

430 PRINT BVOLT 

440 PRINT VDC 

450 FOR I = 1 TO NAX; PRINT AX< I ) : NEXT 
470 PRINT D^; -CLOSE -;FS 

480 HOME : VTAB 15; INPUT "MOTHER RUN^ (1— YES;0— NO) - ; SB): : IF SBX > = 

1 GOTO 120 

490 PRINT -END OF SAMPLING PROGRAM*: END 


QR!GJf\'?iL PA''£ lo 

OF POOR QUALrTY 


9500- 

16 


CLC 


9578- 

A6 

09 

LDX 

409 

9501- 

AO 

02 

LDY 

4402 

957A- 

91 

5A 

STA 

<45A) ,Y 

9503- 

81 

68 

LDA 

(46B) ,Y 

9570- 

06 


INY 


9505- 

65 

68 

ADC 

468 

957D- 

AD 

DO CO 

LDA 

4C0D0 

9507- 

85 

3A 

STA 

43A 

9580- 

8E 

DO CO 

STX 

4C0D0 

9509- 

08 


INY 


9583- 

91 

SA 

STA 

(45A) ,Y 

950A- 

61 

68 

LDA 

(468) ,Y 

9585- 

08 


IKY 


950C- 

65 

6C 

ADC 

46C 

9586- 

DO 

D2 

8NE 

4955A 

950E- 

65 

38 

STA 

438 

9588- 

E6 

56 

INC 

458 

9510- 

AS 

68 

LDA 

46B 

958A- 

06 

05 

DEC 

405 

9512- 

85 

4A 

STA 

44A 

958C- 

DO 

00 

BNE 

4955A 

9514- 

A5 

60 

LDA 

46C 

958E- 

60 


RTS 


9516- 

85 

48 

STA 

448 

958F- 

00 


BRK 


9518- 

A5 

40 

LDA 

440 

9590- 

AO 

08 

LDY 

#408 

951A- 

85 

18 

STA 

418 

9592- 

B1 

68 

LDA 

(468) ,y 

951C- 

AS 

41 

LDA 

441 

9594- 

85 

05 

STA 

405 

951E- 

85 

19 

STA 

419 

9596- 

A9 

08 

LDA 

#408 

9520- 

A5 

06 

LDA 

406 

9598- 

18 


CLC 


9522- 

85 

09 

STA 

409 

9599- 

65 

68 

ADC 

468 

9524- 

20 

40 95 

JSR 

49540 

9598- 

85 

SA 

STA 

45A 

9527- 

18 


CLC 


959D- 

A9 

00 

LDA 

#400 

9528- 

AS 

3A 

LDA 

43A 

959F- 

65 

6C 

ADC 

46C 

952A- 

85 

4A 

STA 

44A 

95A1- 

85 

58 

STA 

458 

952C- 

A5 

36 

LDA 

438 

95A3- 

08 


INY 


952E- 

85 

48 

STA 

448 

95A4- 

81 

68 

LDA 

(468) ,Y 

9530- 

A5 

42 

LDA 

442 

95A6- 

85 

45 

STA 

445 

9532- 

85 

18 

STA 

418 

95A8- 

08 


INY 


9534- 

AS 

44 

LDA 

444 

95A9- 

81 

68 

LDA 

<46B) ,Y 

9536- 

85 

19 

STA 

419 

95AB- 

85 

44 

STA 

444 

9538- 

A5 

07 

LDA 

407 

95AD- 

AO 

00 

LDY 

#400 

953A- 

85 

09 

STA 

409 

95AF- 

81 

5A 

LDA 

<45A> ,Y 

953C- 

20 

40 95 

JSR 

49540 

9581- 

85 

43 

STA 

443 

953F- 

60 


RTS 


9583- 

08 


IKY 


9540- 

AO 

08 

LDY 

8408 

9584- 

Bl 

5A 

LDA 

<45A) ,Y 

9542- 

81 

4A 

LDA 

(44A) ,Y 

9586- 

85 

42 

STA 

442 

9544- 

85 

05 

STA 

405 

9588- 

A5 

42 

LDA 

442 

9546- 

A9 

09 

LDA 

#409 

95BA- 

38 


SEC 


9548- 

18 


CLC 


95BB- 

E5 

44 

SBC 

444 

9549- 

65 

4A 

ADC 

44A 

9580- 

65 

40 

STA 

440 

9548- 

85 

5A 

STA 

45A 

95BF- 

A5 

43 

LDA 

443 

954D- 

A9 

00 

LDA 

#400 

95C1- 

E5 

45 

SBC 

445 

954F- 

65 

48 

ADC 

44B 

95C3- 

85 

41 

STA 

441 

9551- 

85 

58 

STA 

45B 

9505- 

A5 

43 

LDA 

443 

9553- 

AO 

00 

LDY 

#400 

95C7- 

85 

45 

STA 

445 

9555- 

A6 

09 

LDX 

409 

9509- 

A5 

42 

LDA 

442 

9557- 

BE 

DO CO 

STX 

4CODO 

9508- 

85 

44 

STA 

444 

955A- 

AS 

18 

LDA 

418 

95CD- 

88 


DEY 


955C- 

85 

18 

STA 

41B 

95CE- 

AS 

41 

LDA 

441 

955E- 

AS 

19 

LDA 

419 

95D0- 

91 

SA 

STA 

<45A> ,Y 

9560- 

85 

lA 

STA 

41A 

95D2- 

08 


IKY 


9562- 

AS 

lA 

LDA 

41A 

95D3- 

AS 

40 

LDA 

440 

9564- 

FO 

04 

6EQ 

4956A 

95D5- 

91 

SA 

STA 

(45A) ,Y 

9566- 

C6 

lA 

DEC 

41A 

95D7- 

ce 


IKY 


9568- 

DO 

F8 

BNE 

49562 

95D8- 

DO 

D5 

BNE 

495AF 

956A- 

AS 

18 

LDA 

418 

95DA- 

E6 

58 

INC 

456 

956C- 

FO 

07 

BEQ 

49575 

95DC- 

06 

05 

DEC 

405 

956E- 

06 

18 

DEC 

418 

95DE- 

DO 

OF 

BNE 

495AF 

9570- 

06 

lA 

DEC 

41A 

95E0- 

60 


RTS 


9572- 

40 

62 95 

JMP 

49562 

95E1- 

DO 

CO 

BNE 

495A3 

9575- 

AD 

D1 CO 

LDA 

40001 

95E3- 

91 

SA 

STA 

<45A) ,Y 






95E5- 

08 


INY 







95E6- 

DO 

D2 

BNE 

495BA 






95E8- 

£6 

58 

INC 

45B 






95EA- 

06 

05 

DEC 

405 






95EC- 

DO 

00 

BNE 

495BA 






95EE- 

60 


RTS 



B. FIXX.FOR 


FIXX is the FORTRAN program which takes the modified 
data from the APPLE lie data collection program and restores 
it to the full data set by a successive addition process. 

The output data set is used by the following programs AUTOCO 
and FOURIE. 


i 



non onnnrsnno 


Of POOR QUALITY 


PROGRAM FIXX.rOR 


THIS PROGRAh REAlJS IN THE. h(3LiiriE:i« A/D VOLTAGL NUMDERS 
AND RECREATES THE ORIGINAL DATA FILE FOR PROCESSING PV 
EITHER AUTOCO*FOR OR FOURIE*FOR 

DIMENSION NPUT< 20000) 

READ(lf*)NPTS 

WRITEC22f*)NPTS 

READ(lf»)IGAIN 

WRITE(22f*)IGAIN 

READ<lf*)DT 

WRITE<22f»)DT 

READClflOVMKS 

WRITE(22f »)OMKS 

READ(lf»)VDC 

WRITE(22f*)0DC 


THE FIRST VOLTAGE NUMDER IS UNMODIFIED* 

READ(lr20)NPUT(l) 

WRITE(22f20)NPUT(l) 

DO 40 I = 2fNPTS 
J « I-l 


20 

READ( 1 f20fERR=50) IX 
FORMAT(IG) 

40 

NPUT(I) = IX + NPUT(J) 
WRITE(22f20)NPUT( I > 
CONTINUE 

50 

STOP 


END 


C. AUTOCO.FOR 


This FORTRAN program is used on the DEC 1099 computer 
to convert the raw velocity data collected on the APPLE lie 
computer to a more usable form. The output of this program 
is the integral scale and the turbulence intensity. 


nnnno nnnn non *-* »-^nnon nnnnoo nnn nnnnnnn 


ORWrNAL PAGE tS 
OF POOR QUALITY 


PROGRAM AUTOCO.FOR 


THIS PROGRAM lAKEIS THC RE‘-CONST ITUTtH DATA COLLPCTPD 
E«Y THt apple: COMPUTEIR 10 PtRFORM THE! AUl UCORRt LAT ION 
AND CALCULATE THE INTEGRAL SCALE* 

IMPLICIT DOUBLE PRECISION (A-HfO-2) 

D I MENS I ON VELOC (20000 ) t RU ( 500 ) f T ( 50P > 

XMAX IS THE FULL SCALE RANGE OF THE A/D CONVERTER 

XMAX - 4096* 

WRITE(6f 8) 

F0RMAT(/5Xf "NUMBER OF PASSES ? ") 

SET THE MAXIMUM NUMBER OF TIME STEPS TO TAKE* MUST BE 
LESS THAN NPTS/2 

READ(5»*)NPASS 
WRITE(6f 12) 

INPUT THE CALIBRATION COEFFICIENTS USED TO CONVERT 
FROM VOLTAGE TO VELOCITY 

F0RMAT(/5Xf 'INPUT VELOCITY CORRELATION COEF F " S ' r /2X » ' A 

READ(5f*)A 

WRITE(6f 13) 

FORMAT <2Xf 'B = ? ' ) 

READ(5f*>B 

READ(lf*)NPTS 

READ(lf»)IGAIN 


DETERMINE THE FULL SCALE VOLTAGE USED ON THE A/D 

IF(IGAIN*EQ*0*0P*IGAIN.E0*4) FS = 5000, 

IF( IGAIN*EQ* J *0R* IGAIN^EQ.b) fS = 1000. 

IF( IGAIN.E0.2*0R* IGAIN.E0.6) FS = 500. 

IF< IGAIN*E0*3*0R* IGAIN*EQ*7) FS = 100. 

IF A BIPOLAR RANGE WAS SELECT ED r REDUCE THE KESOLUTIUN 
BY HALF* THIS ALSO SETS THE ZERO AT HALE- SCALE* 

IF(IGAIN*GT*3)XMAX = 2048. 

READ(lf*)DT 

READ(lf*)VMKS 

READ(1^*)VDC 


CALCULATE THE WIND TUNNEL MEAN VELOCITY FROM A 
CONSIDERATION OF THE DYNAMIC PRESSURE AND THE MEASURED 
VOLTAGE FROM THE MKS PRESSURE TRANSDUCER* 

A1 = 1*049#VMKS 

B1 = *03001lJ|f(VMKS*»0. 999708) 

C = 531**(A1 - Bl) 

D = 29*6*2*54 
E = D - A1 
F = (C/E)**0*5 
U = 58.582*F 


c 

C IF THE ANEHOMETEF VC ^^OLTAGE IS LESS THAN 3*1 VOLTS 

C THEN THE PROBE WAS IN THE LOU SPEED SECTION* CORRECT 

C THE MEAN VELOCITY AS MEASURED BY THE MKS ACCORDINGLY* 

C 

IF ( VDC*LE*3100* )U = U*C.149 
WRITE(6f 14) 

14 FORMAT (1 HI f 10Xf60(lH») ) 

WRITE<6> 15)U 

15 FORMAT(/24Xf 'MKS VELOCITY = 'rF10*4f' FT/SEC ') 

C 

C READ IN THE A/D VOLTAGE NUMBERS AND CONVERT TO A VOLTAGE 

C 

DO 25 I = IrNPTS 
READ(lf20fERR=30)IX 
20 FORMAT(IG) 

X FLOAT(IX) 

IF(IGAIN*LE*3)G0 TO 24 
C 

C CORRECT FOR A BIPOLAR RANGE 

C 

X * X - XMAX 

24 VELOCd) = X*FS/XMAX 

SUM *= SUM + VELOCd) 

25 CONTINUE 
GO TO 40 

30 NPTS = I 

C 

C CALCULATE THE AVERAGE VOLTAGE* IF NOT ZERO» SUBTRACT 

C FROM EACH VOLTAGE* 

C 

40 VAVG « SUM/FLOAT(NPTS) 

SUM » 0* 

DO 50 I = ItNPTS 
VDIFF = VELOCd) - VAVG 
C 

C ADD EACH ZERO-BASED VOLTAGE TO THE MEASURED DC VOLTAGE 

C AND CONVERT TO VELOCITY USING CORRELATION COEFF 1 C lENTS * 

C 

VOLT = (VDC + VDIFF)/1000* 

VELOCd) = A*(VOLT**B) 

50 SUM = SUM + VELOCd) 

VAVG = SUM/FLOAT(NF‘TS) 

C 

C CALCULATE THE RMS VOLTAGE TO GET THE INTENSITY 

C 

SUM = 0* 

DO 60 I * If NPTS 

SUM * SUM + (VELOCd) - VAVG)**2* 

VELOCd) = VELOCd) - VAVG 
60 CONTINUE 

VRMS = SORT ( SUM/FLOAT ( NPTS ) ) 

TI = VRMS/VAVG*100* 

C 

C PERFORM THE AUTOCORRELATION ON THE FLUCTUATING 

C COMPONENTS OF THE VELOCITY* 

C 

DO 90 I = IfNPASS 
SUM = 0* 

DO 70 J = lfNPTS-d-1) 

SUM * SUM + VELOC(J)*VELOC( J+I-1 ) 

70 PT = FLOAT(J) 
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C 

c 

c 


c 

c 

c 


80 

90 


RV IS THE AUTOCORRELATION COEFFICIENT 

RV(I) *= <SUM/PT)/(VRMS*VRMS) 

T(I> = DT*FLOAT( 1-1) 

END THE PROCEDURE IF RU<I) IS LESS THAN ZERO; 

IF(RV(1) .LT.O. >00 TO 95 
URITE(21»B0> T(l>iRV(l) 

FORMAT <2X»F12,A»SX*F12. 6) 

CONTINUE 
60 TO 98 


95 NPASS - I 
RV(NPASS) = 0. 

WRITE(21 » 80) T(NF'ASS)fRV( NPASS) 

URITE(6f 96)NPASS 

96 F0RMAT(/15Xr 'NEGATIVE CORRELATION AFTER '»I2r' STEPS') 
URITE(6t97)T(NPASS) 

97 FORMAT< 15X. 'LAST STEP WAS 'rF9.3»' MICROSECONDS') 


C 

C THE INTEGRAL SCALE IS FOUND BY USING SIMPSON'S RULE 

C TO INTEGRATE THE AREA UNDER THE RV VS. T CURVE 

C 

98 H = DT*l,0E-06 

SINT = 0. 

DO 100 I = 2»NPASS-1 
100 SINT * SINT + RV(I) 

TSUM = H*SINT + (H/2. )*(RV(1) + RV(NPASS)) 


WRITE(6.105)VAVGfVRMS 

105 F0RMAT(//24Xf 'AVERAGE VELOCITY = '.F12.5.' FT/SEC' f/26X 



C f'HHS VELOCITY * FT/SEC') 

URITE(6f 108>TI 


108 

F0RMAT(27Xf 'TURBULENCE INTENSITY = / 

WRITE(69ll0)TSUM 

i ' ) 

110 

F0RMAT(//29X, 'INTEGRAL = '»1PE12.5»' SEC') 
VSUM = VAVG^TSUM»12. 

WRITE(6f 120)VSUM 


120 

F0RhAT(29Xf 'INTEGRAL = 'fF12»6p' IN*') 
WRITE(69l30) 


130 

F0RhAT(//10Xr60< 1H») r ////> 

STOP 

END 



D. FOURIE.FOR 


This FORTRAN program for use on the DEC 1099 computer 
executes the Fast Fourier Transform of the raw velocity data 
collected by the APPLE lie computer. The output results are 
used with a plotting program to examine the power spectrum 
of the data. The IMSL subroutine FFTRC which actually 
performs the Fast Fourier Transform must be linked to the 
program during execution using the command 

EX FOURIE.FOR, PUB: IMSL. REL/SEARCH 


ORIGINAL PAGE IS 
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c 

C PROGRAM POURIE.FOR 

C 

C 

C THIS PROGRAM PERfORMS THE FASl FOURIER 1 RANGE ORM OE 

C THE RE-CONSTITUTED DATA FROM THE APPLE COMPUTER. AN 

C AVERAGING PROCESS IS USED TO SMOOTH THE RESULTS. 

C 

C THE IMSL SUBROUTINE FFTRC IS USED TO PERFORM THE FAST 

C FOURIER TRANSFORM. 

C 

DIMENSION VEL ( 10000 ) -STORE < 10000 ) t SGK ( 2000 ) » IWK ( 1000 ) f UK ( 1000 ) 
COMPLEX FC(5001) f TEMP 

SMALL = l.OE-06 
PI = 3.1415927 
XMAX = 4096. 

URITE(6>10) 

C 

C THE NUMBER OF POINTS USED TO FIND THE FFT MUST BE A 

C POUER OF 2. HERE 2*»11 = 2048 WAS USED. 

C 

10 FORMAT (///5X. 'NUMBER OF POINTS?') 

READ(5f»)N 


READ(lr#)NPTS 

REAIKlrUOIGAIN 

C 

C DETERMINE THE FULL SCALE VOLTAGE RANGE USED ON THE A/D 

C 

IF(1GAIN.EQ.0.0R.IGAIN.EQ.4) FS = 5000. 
IFdGAIN.EQ.l.OR.IGAIN.EQ.S) FS = 1000. 
IF(IGAIN.E0.2.0R.IGAIN.E0.6) FS = 500. 

IF(IGAIN.E0.3.0R.IGAIN.EQ.7) FS = 100. 

C 

C IF A BIPOLAR RANGE UAS SELECTED. REDUCE THE RESOLUTION 

C BY HALF. THIS ALSO SET THE ZERO AT HALF -SCALE. 

C 

IF(1GA1N.GT.3)XMAX = 2048. 

C 

C DT IS THE TIME BETWEEN EACH SAMPLE. THE MRS VOLTAGE AND 

C THE ANEMOMETER DC VOLTAGE ARE NOT USED HERE. 

C 

READ(1.*)DT 

READ(1»*)VMKS 

READ(1.*)VDC 

C 

C READ IN THE A/D VOLTAGE NUMBERS AND CONVERT TO A VOLTAGE 

C 

SUM = 0. 

DO 25 I = 1 . NPTS 
READ(1.20.ERR=30)IX 
20 FORMAT <1G) 

X = FLOAT(IX) 

IF(I6AIN.LE.3)G0 TO 26 
X = X - XMAX 

26 STORE(I) = X*FS/XMAX 

SUM = SUM + STORE(l) 

25 CON’-INUE 

GO TO 40 
30 NPTS = I 


IF NOT SUKTKACT 


C 

C CALCULATE THE AVERAGE VOLTAGE* 

C FROM EACH VOLTAGE* 

C 

40 VAVG == SUM/FLOAT (NF TS) 

DO 50 I = IrNPTS 
STORE(I) = STORE(I) - VAVG 
50 CONTINUE 

C 

C THE RECORDS USED FOR THE FFT OVERLAP BY 1/2 THE NUMBER 

C OF SAMPLES IN THE RECORD. THIS GIVES A BETTER AVERAGE. 

C 

NOVER := N/2 
NTIMES = 2*(NPTS/N) 

NSTART * 0 

DO 60 I = If 2000 
C 

C ZERO THE ARRAY OF AVERAGED FOURIER COEFFICIENTS 

C 

60 SGK(I) = 0. 

DO 140 NCOUNT = IfNTIMES 
C 

C MAKE SURE THE RECORD HAS A ZERO AVERAGE. 

C 

SUM = 0. 

DO 70 I = IfN 

VEL<I) = STORECNSTART +1) 

70 SUM = SUM + VEL(I) 

VAVG == SUM/FLOAT (N) 

DO 80 I = IfN 

80 VEL(l) = VEL(I) - VAVG 

C 

C USE A HAMMING WINHOU ON THE FIRST AND LAST 10/1 OF THE 

C RECORD TO REDUCE THE APPEARANCE OF SIDELOBES CAUSED 

C BY A FINITE RECORD LENGTH. 

C 

M = N/10 
IFRONT = M 
IBACK = N-M 

DO 100 I = IfN 
IFd.GE. IFRONT) GO TO 90 
J = I“1 

VAR = C0S(PI»FL0AT(J)/FL0ATCM-1) > 

VEL(I) = (VEL(I)/2. )»(1. - VAR) 

GO TO 100 

90 IFd.LE. IBACK) GO TO 100 

J = N-I 

VAR = COS(PI»FLOAT( J)/FL0AT(M-1 ) ) 

VEL(I) = <VEL(I)/2.)*(1. - VAR) 

100 CONTINUE 
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c 

C PERFORM THE FAST FOURIER TRANSFORM ON THE RtCORl' 

C 

CALL FFTRC(VELtN»FC»IWKfUK) 

C 

C AVERAGE EACH DATA POINT IN THE RECORD WITH THE DATA 

C POINTS ON EACH SIDE* 

C 

DO 110 I = 2f (N/2) 

TEMP ' 0*25*FC(I-1) + 0*5»FC(I) + 0.25*FC<I+1) 

110 FC(I) - TEMP/0*875 

C 

C NORMALIZE THE COMPLEX FOURIER COEFFICIENTS BY COMPUTING 

C THE ABSOLUTE VALUE AND DIVIDING BY THE NUMBER OF SAMPLES 

C IN THE RECORD* COMPUTE THE AVERAGE FOR EACH RESULTING 

C FREQUENCY 

C 

FREQ = 0* 

GDC = 10**AL0G10( < (CABS(FC( 1 ) )**2. )/SMALL )/FL0AT<N/2) ) 
SGK(l) = SGK(l) + GDC 

J = N/2 

NSTEP = J/1000 

IF(NSTEP*LE*1> NSTEP = 1 

DF = 1*/(2*»FL0AT( J)*(5MALL*DT) ) 

NGK = 1 

DO 130 I = If (J/NSTEP) 

IGK = I + 1 
FSUM 0* 

GSUM * 0* 

DO 120 L * If NSTEP 
FREQ = FREQ + DF 
NGK = NGK 4 1 

6K = 10**AL0G10( < (CABS(FC(NGK> )**2* )/SMALL)/FLOAT( J) ) 
FSUM = FSUM + FREQ 
GSUM = GSUM 4- 6K 
120 CONTINUE 

GAVG = GSUM/FLOAT (NSTEP) 

SGK(IGK) = SGK(IGK) 4 GAVG 
130 CONTINUE 

NSTART * NSTART + NOUER 

lAO CONTINUE 

C 

C OUTPUT THE AVERAGED FREQUENCIES AND NORMALIZED 

C SPECTRUM AT THAT FREQUENCY IN DB* 

C 

FREQ = 0. 

DO 150 I = lf<J/NSTEP) 

GAVG = SGKd )/FLOAT(NTIMES> 

WRITE(22f *)FREQfGAVG 

FREQ * FREQ 4- DDK ( FLOAT < NSTEP ) ) 

150 CONTINUE 


STOP 

END 


E. REDUCE. FOR 


REDUCE. FOR for the DEC 1099 computer takes the data 
produced by AUTOCO.FOR for one velocity and all grid posi- 
tions and calculates the correlations for the integral scale 
and turbulence intensity. It also calculates the dissipation 
length parameter by using the turbulence decay law proposed 
by Castro. A nonlinear optimization method is used to find 
the virtual origin which gives the best power-function fit 
to the integral scale results. 



non 


PROGRAM REDUCE ♦FOR 


C 
C 
C 

C THIS PROGRAM TAKES THE GRID POSITION^ INTEGRAL SCALEr 

C AND TURBULENCE INTENSITY DATA AT ONE MEAN VELOCITY AND 

C PERFORMS A CORRELATION FOR THE INTEGRAL SCALE AND 

C THE DISSIPATION LENGTH PARAMETER. 

C 

IMPLICIT DOUBLE PRECISION (A-HtOtZ) 

DOUBLE PRECISION NAME 

DIMENSION XGRID(20) r TSCALE ( 20 ) » XSCALE < 20 ) r TPRIME < 20 ) r XM ( 20 ) 
DIMENSION X<20) f Y(20) f XD(20) f YD(20) 

COMMON /X/X/XSCALE/XSCALE/AX/AX/BX/BX/RX/RX/NPTS/NPTS 

IlOA = l.OE-04 
Vise * 15.69E-06 
NAME * 'SMALL' 

READ (lr*)NPTS 
READ <lf»)XMESH»B 
READ <1»*)XPR0BE 
READ <lf»>SPEED 
IF(SPEED.LT.20. )NAME = 'LARGE' 

DO 10 I = 1»NPTS 

REALUlf*)XGRID(I)fTSCALE(I)fTPRIME(I) 

10 CONTINUE 


SORT THE DATA IN ORDER BY GRID POSITION 

XN * FLOAT(NPTS) 

ND = 2*»(ALOG(XN)/ALOG(2. ) ) - 1 
20 DO 60 I * Ir(NPTS-ND) 

IF(XGRID(I) .LE.XGRID(I4ND) ) 60 TO 60 
TX * XGRID(I+ND) 

TLT * TSCALE(I+ND) 

TTURB = TPRIMECI+ND) 

XGRID<I+ND) = XGRID(I) 

TSCALE(I+ND) = TSCALE(I) 

TPRIME<I+ND) = TPRIML(I) 

IF(I.GT.ND) GO TO 30 
XGRID(I) = TX 
TSCALE(I) = TLT 

TPRIME(I) = TTURB 

GO TO 60 

30 DO 40 J = (I-ND) ^ (-NIO 

IF(TX.GE.XGRID< J) > GO TO 50 
XGRID(J+ND) - XGRID(J) 

TSCALECJ+ND) « TSCALE(J) 

40 TPRIMEC J+ND) = TPRIME(J) 

50 XGRID<J+ND) * TX 

TSCALE( J+ND) = TLT 
TPRIME(J+ND) = TTURB 
60 CONTINUE 

ND = ND/2 

IF(ND.GT.0)G0 TO 20 
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CALCULATE THE INTEGRAL SCALE FROH THE TIME SCALE HATA 
CALCULATE THE TURBULENT ENERGY DECAY < U**2 ) / ( U ' **2 ) 

DO 70 I = l.NPTS 

X(I> * XGRID(I) - XPROBE 

XM(I) = X<I)/XMESH 

XSCALE(I) * <TSCALE(I)*T104*SPEED)»100. 

Y(I) = (TPRIME<I)/100.)**2. 

YD< I ) = 1 ./Y( I ) 

CONTINUE 


FIND THE DISSIPATION VIRTUAL ORIGIN BY A LINEAR FIT OF 
THE TURBULENCE DECAY DATA 

CALL LFUNC<XrYDrAL»BL»RLfNPTS) 

XOD = -AL/BL 

DO 71 I = l.NPTS 

XD(I) = (X(I) - X0D)/XMESH 


DO A POWER FUNCTION FIT TO THE TURBULENCE DECAY DATA 
USING THE VIRTUAL ORIGIN TO GET THE COEFFICIENTS USED 
TO CALCULATE THE DISSIPATION PARAMETER. 

CALL PFUNC < XD . Y . AD . BD . RD . NPTS ) 

CALL PFUNCCX.TPRIMEiAT.BT.RTf NPTS) 


USE A UNIVARIATE LINE SEARCH ALGORITHM TO FIND THE BEST 
FIT TO THE INTEGRAL SCALE DATA USING A VIRTUAL ORIGIN 

NVAR = 1 
XO = 0. 

CALL FUNC(XO.FXO.NVAR) 

CALL GRADIE<XO.FXO. FART. NVAR) 

SMAG = SQRT(PART»*2. ) 

TEST = -PART/SMAG 

CALL GOLDENCXO. TEST. XORG. NVAR) 

WRITE<6.72)XMESH.NAME 
72 F0RMAT(1H1.20X. 'DATA FOR '.F5.2.' cm GRID IN '.A5 

C . ' TEST SECTION' ) 

WRITE (6. 74) SPEED. XPROBE 

74 F0RMAT(/33X. 'VELOCITY = '.F5.2.' m/s '. /33X .' PROBE AT X = ' 

C F6.2. ' cm' .//) 

WRITE (6. 80) .XORG.RX.AX.BX 

80 FORMATC lOX. ' INTEGRAL SCALE VIRTUAL ORIGIN AT Xo = '.F7.2.' 

C /20X.'R = ' .F9.6./20X. 'A = ' . 1PE12 . 4 . /20X . ' B = ' . 1PE12 . 4 . // ) 
WRITE(6.90)X0D.RD.AB.BD 

90 FORMATdOX. 'DISSIPATION VIRTUAL ORIGIN AT Xo = '.F7.2,' cm' 

C /20X.'R = ' .F9.6./20X. 'A = ' . IPEl 2 . 4 . /20X . ' B = ' . IPEl 2 . 4 . // ) 
WRITE(6.94)RT.AT.BT 

94 FORMAT( lOX. 'CORRELATION FOR TURBULENCE INTENSITY IS'./20X. 

C 'R = ' .F9.6./20X. 'A = ' . 1 PE12 . 4 . /20X . ' B = ' . 1PE12 . 4 . /// ) 


? 


chi ' f 


WRITE(6f 100) 

100 FORMA! (//T8* 'GRID ' rT17r ' X ' t T22f 'X-Xoi ' fT29f ' X-Xod ' f T36 ^ ' TURD 

C T45r 'TIME' #T54r 'INT' fJ63f 'DIGS' fT73r 'Lx' > 

WRITE<6tllO) 

110 FORMAT (T6r 'POSITION ' »T16 f ' ' »T22r ' ' tT29f '* ' r 

C T35f 'INTENS' fT44f' SCALE ' fT53f 'SCALE' »T62f 'SCALE' fT72f ' ' ) 

WRITE(6f120) 

120 FORMAT (TBf ' (cit.) ' t T17» 'M ' t T24 1 'M ' t T31 f 'M' tT37f ' i%) ' tT45f ' (uS) 

C T54f ' (cii.) ' rT63f ' (cm) ' rT73f 'Ld' »/T4f 74( 1H=) ) 


RB « SPEED*B/<100.*0ISC) 

DO 140 I *= IrNPTS 
XG « XGRID(I) 

TP = TPRIME(I) 

C 

C CALCULATE THE DISSIPATION SCALE USING THE VIRTUAL ORIGIN 

C AND THE POWER-FUNCTION COEFFICIENTS* 

C 

DS = (AD**0*5)*XMESH»(XD(I)»t(l. + BD/2 • > ) / < -BD ) 

XTD = XSCALE(I)/DS 
TM « TSCALE(I)*100* 

XOBR = (X<I) - XORG)/(B*RB) 

YBR - XSCALE(I)/(B> 

WRITE<6f 130)XGFXM(I)fX0MfXD(I) fTPfTMfXSCALEC I) f DSfXTD 
130 F0RMAT(T7fFA.2tT15FF5*2pT21fF6*2fT28fF6*2fT36fF5*3f T43iF7.2t 

C T53fF6*4FT62rF6.4tT71fF6*4) 

HRITE(30r*)X( I ) fTPRIMEC I ) 

WRITEC31 f *)X( I ) fXSCALEC I ) 

WRITE(32»*)X< I) fDS 
WRITE(9f<t)X0BRfYBR 
140 CONTINUE 

WRITE<6f 150) 

150 FORMAT (IHl) 

STOP 

END 


C 

C GOLDEN IS THE UNIVARIATE LINE SEARCH ALGORITHM WHICH 

C USES THE GOLDEN SECTION BRACKETING TECHNIQUE 

C 

SUBROUTINE GOLDEN ( X f S r XSTAR » N ) 

IMPLICIT DOUBLE PRECISION ( A-Hf O-Z ) 

FI = 0*381966 
F2 = 0*618034 
STOL = 0*001 

CALL FUNC(XfFXfN) 

FXL = FX 
XL = X 

DO 10 I = 1f800 

XH = XL + S 

CALL FUNC(XHfFXHfN) 

IF(FXH*GT*FXL) GO TO 20 
XL = XH 
FXL FXH 
10 CONTINUE 



20 IF( I *E0. 1 ) GO TO 30 

XL = XH - 2**S 
CALL FUNC(XLrFXL»N') 

30 DELTA * XH " XL 

DO 60 IB = 1»300 

DELI = F1*DELTA 
DEL2 = F2*DELTA 
XI = XL + DELI 
X2 = XL + DEL2 
CALL FUNC(X1 fFXl f N) 

CALL FUNC(X2f FX2»N> 
IF<FX2.GE*FX1) GO TO 40 
DELTA » XH - XI 
XL ie XI 
GO TO 50 

40 DELTA = X2 - XL 

XH = X2 

50 IF( ABS(DELTA) ♦LE*STOL) GO TO 70 

60 CONTINUE 

70 XSTAR * X2 

RETURN 
END 


C 

C THIS SUBROUTINE IS USED BY GOLDEN AS THE OBJECTIVE FUNCTION 

C THE TASK IS TO MAXIMIZE THE POWER -FUNCl ION CORRELATION 

C COEFFICIENT RX BY CHOICE OF THE VIRTUAL ORIGIN XO * 

C 

SUBROUTINE FUNC ( XO f FX f NVAR ) 

IMPLICIT DOUBLE PRECISION (A-HfO~Z) 

DIMENSION Y(20) fTRY(20) f X < 20 ) f XSCAL E ( 20 ) 

COMMON /X/X/XSCALE/XSCALE/AX/AX/BX/BX/RX/RX/NPTS/NF TS 

DO 10 I = IfNPTS 
TRY<I) = X(I ) - XO 
Yd) XSCALE(I) 

10 CONTINUE 

CALL PFUNCCTRYFYfAXfBXFRXFNPTS) 

FX = -<RX) 

RETURN 

END 


C 

C PFUNC IS THE POWER-FUNCTION LEAST-SQUARES ANALYSIS FUNCTION 

C 

SUBROUTINE PFUNC ( Xf Y f A f B fR f N ) 

IMPLICIT DOUBLE PRECISION (A-HfO-Z) 

DIMENSION X<20) f Y(20) 



TSMALL «= 0*000001 
XI = 0* 

Y1 = 0* 

SI = 0* 

T1 = 0* 

W1 = 0. 

DO 10 I * 1»N 
IF(X(I) .LE*0* )G0 TO 2 
XL * DLOG(Xd) ) 

GO TO 4 

XL * DLOG( TSMALL) 

IF(Y(I).LE.O.:)GO TO 6 
YL = DLOG(Y<I)) 

GO TO 8 

YL = DLOG<TSMALL) 

XI * XI + XL 
Y1 « Y1 + YL 
SI ■= SI + XL*YL 
T1 » T1 + XLtXL 
W1 = W1 + YL»YL 

XN = FLOAT(N) 

D1 = XN*T1 XltXl 
AL = (Y1*T1 - X1*S1)/D1 
B = (XN3ICS1 - X1*Y1)/D1 

R = (AL»Y1 + B*S1 - Y1*Y1/XN)/(W1 - Y1»Y1/XN> 
A = DEXF(AL) 

RETURN 

END 


LFUNC IS A LINEAR LEAS7 -SQUARES ANALYSIS ROUTINE 


SUBROUTINE LFUNC C X f Y fA. B f R f N ) 
IMPLICIT DOUBLE PRECISION <A-H.O-Z> 
DIMENSION X(20)»Y(20> 


XI 


0. 


Y1 


0* 


SI 


0. 


T1 

= 

0* 


U1 

= 

0. 


DO 

10 I 

= 

XI 


XI 

+ 

Y1 

= 

Y1 

+ 

SI 

= 

SI 

+ 

T1 

=5 

T1 

+ 

W1 

= 

kll 

+ 


1 rN 
X(I ) 

Y(I) 

X(I)*Y(I) 

X<I)»X(I) 

Y<I)*Y(I) 


XN = FLOAT(N) 

Dl = XN*T1 - X1*X1 
A = <Y1*T1 - X1»S1 )/Dl 
B = (XN*S1 - X1*Y1)/U1 

R = (A*Y1 + B*S1 - YUY1/XN)/(W1 - Y1*Y1/XN> 


RETURN 

END 
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GRAHIE IS USED TO FIND THE DIRECTION IN WHICH TO BEGIN 
SEARCHING TO ACHIEVE THE BEST FIT FOR THE INTEGRAL SCALE 
CORRELATION* 

SUBROUTINE GRADIE(XfFX»PART»N) 

IMPLICIT DOUBLE PRECISION (A-HrO-Z) 

DX *= l*0E-04 

XI = X + DX 

CALL FUNC(XlfFXlrN) 

part * (FXl - FX)/nX 

RETURN 

END 



Test Section 


(Re/L) X 10-^ 



* 

1.02 X 1.02 m 

11.2 m/s 

71.3 m-1 


(40 X 40 in.) 

36.7 fps 

1.81 in.-l 

• 

0.411 X 0.411 m 

74.7 m/s 

462 m”l 


(16.2 X 16.2 in.) 

245 fps 

11.7 in.-l 


Table !♦- VUES wind tunnel flow data* Assumed room 

air at typical conditions of 1.002 x 10^ N/m^ 
(29.60 in.Hg) and 295 K (72°F). No model 
installed in tunnel. 


Grid No./ 


Grid M cm 

Grid X cm 

U 00 m/s 

T’„% 

1^ 00 cm 

00 cm 

#1/10.16 

265.7 

25.6 

0.45 

0.82 

0.65 


234.5 

25.6 

0.53 

0.70 

0.53 


265.7 

50.6 

0.57 

0.98 

0.81 


265.7 

71.3 

0.63 

1.09 

0.86 


234.5 

50.6 

0.66 

0.86 

0.67 


215.5 

25.6 

0.70 

0.64 

0.44 


234.5 

71.3 

0.74 

1.06 

0.74 


215.5 

71.3 

0.85 

0.94 

0.65 


194.0 

25.6 

0.87 

0.64 

0.33 


215.5 

50.6 

0.88 

0.81 

0.58 


194.0 

50.6 

0.99 

0.83 

0.46 


173.2 

25.6 

1.01 

0.57 

0.19 


194.0 

71.3 

1.02 

0.92 

0.54 


173.2 

50.6 

1.16 

0.72 

0.31 


173.2 

71.3 

1.18 

0.86 

0.40 

#2/5.08 

265.7 

25.6 

0.22 

1.43 

0.34 


234.5 

25.6 

0.22 

1.09 

0.31 


265.7 

50.6 

0.26 

1.51 

0.38 


194.0 

25.6 

0.28 

0.84 

0.26 


265.7 

71.3 

0.31 

1.38 

0.56 


194.0 

50.6 

0.1 

0.78 

0.30 


215.5 

71.3 

0.32 

0.81 

0.49 


173.2 

25.6 

0.32 

0.73 

0.23 


194.0 

25.6 

0.37 

0.71 

0.46 


173.2 

50.6 

0.38 

0.70 

0.27 


153.2 

25.6 

0.38 

0.63 

0.20 


173.2 

71.3 

0.40 

0.53 

0.43 


153.2 

50.6 

0.46 

0.45 

0.24 


153.2 

71.3 

0.48 

0.43 

0.39 

#3/5.08 

265.7 

26.4 

0.45 

0.79 



265.7 

47.0 

0.46 

0.59 



265.7 

64.5 

0.50 

0.45 



215.5 

26.5 

0.60 

0.61 



215.5 

51.5 

0.66 

0.58 



215.5 

62.8 

0.74 

0.68 



153.2 

26.4 

1.08 

0.53 



153.2 

50.3 

1.19 

0.58 



153.2 

64.6 

1.17 

0.64 


#4/2.54 

58.4 

25.0 

2.59 

0.64 



58.4 

49.8 

2.73 

0.66 



58.4 

55.8 

2.73 

0.71 


#5/5.08 

58.4 

52.5 

3.16 

0.77 


No Grid 


51.6 

0.16 

1.35 



— 

65.1 

0.34 

1.18 



Table 2.- Summary of free-stream turbulence characteristics. 


*.995 cm 


o 

a 

B.L. Tnp 

u jn/s 

T J% 


Exp . * 

xd 

at 

x=22.1 

at 

x=29.7 

1 

yes 

50.5 

0.16 

1.35 

29.6 

0.46 

0.93 

3 

yes 

ti 

•I 

II 

29.7 

— 

— 

6 

yes 

It 

II 

•• 

29.1 

.60 

— 

9 

yes 

II 

M 

It 

27.9 

— 

— 

1 

yes 

64.2 

0.34 

1.18 

30.2 

.44 

.88 

1 

yes 

50.5 

0.46 

0.59 

29 J 

.46 

.93 

1 

yes 

64.2 

0.50 

0.45 

30.4 

— 

— 

1 

yes 

50.5 

1.19 

0.58 

29.7 

.46 

.93 

1 

yes 

64.2 

1.17 

0.64 

30.2 

— 

— 

1 

yes 

00 

o^ 

2.73 

0.66 

29.7 

.46 

.93 

1 

yes 

51.1 

3.16 

0.77 

30.5 

.46 

.93 

1 

no 

50.5 

0.16 

1.35 

30.1 



1 

no 

64.2 

0.34 

1.18 

30.2 




* - 0.25 cm est. 


Table 3.- Detachment locations shown by liquid film and 
corresponding flow conditions. 


' 5.995 


cm 


T * 7 

i 00 /o 


0.67 


0.56 


0.46 29.7 ± 0.25 

1.19 

2.73 

0.31 

0.38 

2.85 


Table 4.- Effect of free- stream turbulence 

intensity upon detachment location 
with constant turbulence Integral scale 


T«i^ lx<*/ '5.995 c*" 


0.30 


1.64 29.7 ± 0.25 


ft 


1.32 

1.23 

0.63 


Table 5.- Effect of free-stream Integral scale 
upon detachment location with 
constant turbulence Intensity. 


X(j cm 


Run 

a® 

U otn/s 

TU 

Ix * cm 

Exp.* 

Strat. 1 

Strat. 2 

396B 

1 

50.5 

0.16 

1.35 

29.6 

28.3 

29.9 

404A 

3 

It 

ti 

II 

29.7 

27.1 

27.7 

404B 

9 

fi 

II 

II 

27.9 

27.9 

27.9 

374A 

1 

It 

1.19 

0.58 

29.7 

27.4 

30.1 

374B 

M 

64.2 

1.17 

0.64 

30.2 

27.7 

29.9 

363 

■t 

49.8 

2.73 

0.66 

29.7 

28.5 

28.8 

403 

• 1 

51.1 

3.16 

0.77 

30.5 

28.4 

30.0 


♦Uncertainty i 0.25 cm 


Table 6.- Comparison of experimental and 
predicted detachment locations 




Xd cm 


Run 

Exp. 

strat. 2 

Townsend 

396B 

29.6 

29.9 

30.8 

404A 

29.7 

27.7 

28.7 

374A 

29.7 

30.1 

30.8 

363 

29.7 

28.8 

29.26 

403 

30.5 

30.0 

30.7 


Table 7 


Comparison of two prediction methods 











Figure 5 


Flow probe support structure 
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Figure 6. Airfoil and mounting system 
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Figure 7. Fiber-optic surface-sensing device and a hot-wire 
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1.0 2.0 3.0 4.0 5.0 

DISTANCEC .001 IN) X 10’^ 

Variation of hot-wire signal with distance from 
surface of model at U_ = 0 
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Figure 10 


Sketch of airfoil cross section with spoiler 




Grid No* 
1 
2 

3 

4 

5 

*1.87 CM 


M CM 

h CM 

Grid Element 
Cross Section 

10.16 

1.24 

rectangular* 

5.08 

0.32 

ROUND 

5.08 

1.24 

rectangular* 

2.54 

0.64 

ROUND 

5.08 

0.64 

ROUND 


IN STREAMWISE DIRECTION 


Figure 11. Dimensions of turbulence-producing grids 
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Figure 13. Effect of wind tunnel contraction on integral 
scales 



COmRACTION RATIO 

Figure 14. Effect of wind tunnel contraction on relative 
turbulence intensity 
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Figure 15. Cross-stream distribution of with minimum X/M 
grid #5 













(c) 





-PRIME sg AT Y=0. 15 MM 



X CM 

Figure 17. Transition locations on the airfoil with no 
boundary layer trips and no grids 
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Figure 18. Liquid-film data on the separation zone; 
views 


overhead 
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Figure 18 cont'd, 
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Figure 20. 
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□ PRESTON TUBE DATA 

■ INFERRED BY FITTING BOUNDARY LAYER PROFILE 
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Figure 23. 


The effect of turbulence intensity upon skin 
friction 



Figxire 24. Boimdary layer velocity profiles for varying free 
strecim turbulence: x = 22.1 cm 
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Figure 25. Boundary layer turbulence intensity profiles for 
varying free-stream turbulence: x = 22.1 cm 



-CM 


(M 



IX- CM 


Figure 26. Boundary layer integral scale profiles for varying 
free-stream turbulence: x = 22.1 cm 
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Figure 27. Boundary layer velocity profiles for varying free 
stream turbulence; x = 29.7 cm 
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SEE FIG. 27 FOR SYMBOL CODE 
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Boundary layer turbulence intensity profiles for 
varying free-streain turbulence: x = 29.7 cm 


Figure 28 
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